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ABSTRACT 

An extraordinary number of Hubble constant measurements challenges physicists with selec¬ 
tion of the best numerical value. The standard U.S. Nuclear Data Program (USNDP) codes and 
procedures have been applied to resolve this issue. The nuclear data approach has produced the 
most probable or recommended Hubble constant value of 66.2(77) (km/sec)/Mpc. This recom¬ 
mended value is based on the last 25 years of experimental research and includes contributions 
from different types of measurements. The present result implies (14.78±1.72) x 10® years as a 
rough estimate for the age of the Universe. The complete list of recommended results is given 
and possible implications are discussed. 

Subject headings: cosmology: observations — methods: data analysis 


1. Introduction 

The Hubble constant and its present-day nu¬ 
merical value play an important role in mod¬ 
ern astrophysics a nd cosmology (Boyd. 2008t 


their policies are well d escribed in lit erature 


Dolgov et al. Ill98^ . For many years, researchers 
have been improvin g the accuracy of the constant 
(jLivio. Riess Il201.lh . Fig. [T] shows the time evo¬ 
lution of Hubble research in the last 100 years, 
and t he original effort can be traced as early as 


1916 (jHuchra et al. 120151) . The large number 


of measurements creates a certain degree of con¬ 
fusion about Hubble constant numerical value, 
and s cientists often rely on recent ly-published re¬ 
sults ( Particle Data Group 112013 ). The precision 
of Hubble constant measurements has improved 
enormously over the years; however, it is not al¬ 
ways prudent to reject older results in favor of the 
latest findings. Consequently, it makes perfect 
sense to analyze all available results, evaluate the 
data, and extract the recommended value. 

2. Hubble Constant Evaluation 

The volume of Hubble constant measure¬ 
ments far exce eds other experim ental quanti¬ 
ties in physics ( Pritvchenko 20151) . Previously, 
similar situations have been encountered in nu¬ 
clear and particle physics and resolved with 
data evaluations. Nuclear data evaluations and 


( Pritvchenko et al. 1 12013: IeI^DF 1 1201511 . Fre 


quently, the evaluations are completely based on 
or adjusted to available experimental data. These 
adjustments and specialized mathematical statis¬ 
tics techniques can be applied for nuclear, particle, 
or any other data sets. 

In this work, 1 would follow standard nuclear 
data evaluation procedures to deduce the recom¬ 
mended value. Current evaluation input data are 
mostly based on the NASA/HST Key Project 
on the Extragalactic Distance Scale compilation 
( Huchra et al. 20151) and recent results. A vi¬ 
sual inspection of historical Hubble Constant mea¬ 
surements, as shown in Fig. [I] is instrumental 
in the data analysis. It suggests that one may 
safely reject all measurements prior to 1970. It 
is common knowledge that Hubble constant mea¬ 
surements heavily rely on the accuracy of astro¬ 
nomic distance determination. Older results, such 
as those reported by A. Sandage and G. de Van- 
couleurs, suffered f rom inaccurate measurements 


(|Livio. Riess I l201.ll ). Therefore, the rejection of 


all results prior to 1990 could provide a compli¬ 
mentary benchmark value of the Hubble constant. 

In the present data analysis, the experimental 
data have been separated into two groups, with 
1970 and 1990 time cuts, and further reduced us- 
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ing the following policies: 

• Rejection of repeated results (multiple pub¬ 
lication of the same result) 

• Rejection of model-dependent results (i.e. 
Cosmic Microwave Background (CMB) fits) 

• Rejection of pot ential outliers using Chau- 
venet’s criterion (jBirch. Singh 1120141 1 


Common data evaluation practices indicate 
that recommended value should be based on 
a large statistical sample that includes differ¬ 
ent types of measurements. The 1970 and 1990 
redacted data sets of ~334 and ~266 data points, 
respectively, provide such samples. These large 
samples create the possibility of deducing Hub¬ 
ble constant value for each method of observa¬ 
tion besides the combined value that is based on 
all measurements. The current data collections 
were further subdivided using a NASA/HST Key 
Project on the Extragalactic Dista nce Scale clas¬ 
sifica tion of experimental methods ( Huchra et al. I 

mi): 


S = Type la supernovae (SNIa) 

2 = Type II supernovae (SNII) 

L = Lens 

r = Red Giants 

B = Baryonic Tully-Fisher 

R = Inverse Tully-Fisher (ITF,RTF) 

H = Infrared Tully-Fisher (or IRTF) 

F = Fluctuations 

A = Global Summary 

Z = Sunyaev-Zeldovich 

T = Tully-Fisher 

O = Other 


These experimental data sets have been pro¬ 
cessed with th e latest version of t he visual aver¬ 
aging library (IBirch. Singh 1201411 . The library 


(FVM) statistical methods to calculate averages of 
experimental data with uncertainties. The exper¬ 
imental data sets were processed, and evaluated 
values with reduced x^<2 were typically accepted 
as reasonable data fits. The current evaluation in¬ 
corporates statistical methods based on the inverse 
squared value of the quoted uncertainties, a proce¬ 
dure that is consistent with the general method¬ 
ology use d in treatment of data for the FNSDF 
database (IfNSD F 112015) an d Particle Data Group 


(jParticle Data Group 11201411 . 


3. Results and Discussion 

Two sets of recommended values are displayed 
in Fig. [21 and the combined numerical values are 
shown in Table |T] The Hubble constant com¬ 
bined central values extracted by means of differ¬ 
ent mathematical techniques are in good agree¬ 
ment, while uncertainties need further discussion. 
Visual inspection of the numerical values shown 
in the Fig. |2] indicates the Best Representation 
mathematical procedure (depicted as “BR”) pro¬ 
vides a best fit for different experimental meth¬ 
ods with reasonable uncertainties. Other proce¬ 
dures such as the Bootstrap technique (depicted as 
“Bs”) result in rather small uncertainties. These 
small uncertainties are due to specifics of data 
analysis procedures such as assignment of different 
statistical weights to results with different uncer¬ 
tainties. In light of this disclosure, it is prudent 
to select the Best Representation results as the 
recommended value. The current results are con¬ 
sistent with the recent Particle Dat a Group pub¬ 
lication ( Particle Data Group 20I4h . 


program includes limitation of relative statisti¬ 
cal weight (LWM), normalized residual (NRM), 
Rajeval technique (RT), and the Expected Value 


Finally, two different time cuts of 1970 and 1990 
for Hubble’s data have yielded two recommended 
values of 66.2(89) and 66.2(77) (km/sec)/Mpc, re¬ 
spectively. The agreement between these values 
partially reflects the fact that the majority of the 
Hubble’s constant measurements has been per¬ 
formed in the last 25 years, and a small number 
of potential outliers has been rejected. More ac¬ 
curate recent observations imply a preference for 
the 1990 time cut value of 66.2(77) (km/sec)/Mpc. 
The last result is consistent with the Hubble Space 
Telescope and Wide Field Gamera 3 ( Riess II 2 OIII I 
and the model -dependent P l anck’s Mission and 
WMAP values ( Bucher II 2 OI 3 I : [Bennett |[2013 1. In¬ 
clusion of the globular cluster luminosity function. 
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Dn-Sig/Fund Plane, Novae and planetary nebula 
luminosity function data would slightly change 
1970 and 1990 recommended values to 66.9(90) 
and 66.9(78), respectively. These less precise mea¬ 
surements do not affect the recommended values 
severely because of a Chauvenet’s criterion analy¬ 
sis. 


In recent years, an effort has been made to 
calculate the Hubb le constant median statistics 
( Chen. Ratra ][20nl)- The median statistics ap¬ 
proach differs substantially from the nuclear or 
particle data evaluation procedures. At the same 
time, it provides a compl ementary value of 6 8T5.5 
(or ±1) (km/sec)/Mpc ( Chen. Ratra 2011 1 that 
can be compared with the present results. 

The knowledge of Hubble constant value has 
multiple implications in science. As an exam¬ 
ple, a rough estimate of the age of the Uni- 
verse can be dedu ced using the standard methodic 
( Wikipedia I I 2 OI 5 II . The adopted value of 66.2(77) 
(km/sec)/Mpc implies (14.78±1.72) x 10® years es¬ 
timated value for the age of Universe. The last re¬ 
sult is consistent with the rece ntly published value 
of (13.798±0.037)xl0® years (|P1a.Tick Il20l 3l . 


4. Conclusions 

The analysis of Hubble constant measurements 
has been performed using standard USNDP codes 
and procedures. An evaluated data set of most 
probable values of Hubble constant has been de¬ 
duced and shown in the Table [TJ These values are 
consistent with other available results. An accu¬ 
rate constant value is instrumental for many po¬ 
tential applications. The recommended value of 
the constant is completely based on experimental 
measurements, and further, more precise observa¬ 
tions, would lead to more accurate determination 
of it. 
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Fig. 1.— Historical evolution of Hubble constant measurements. 
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Fig. 2.— Evaluated values of Hubble constant using 1970 and 1990 time cuts for experimental data. These 
plots include the evaluated or recommended values for combined observations, using the Best Representation 
(BR) and Bootstrap (Bs) data analysis techniques. Globular cluster luminosity function (G), Dn-Sig/Fund 
Plane (D), Novae (N), Planetary nebula luminosity function (P) results are also included while GMB Fit 
(G) value is shown for comparative purposes. 
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Table 1: Results of the Hubble Constant evaluation for all observations using 1970 and 1990 time cuts. 


Method 

Time Cut: 1970 
(km/sec)/Mpc 

Time Cut: 1990 
(km/sec)/Mpc 

Unweighted Average 

67.56(73) 

65.92(65) 

Weighted Average 

63.68(58) 

63.87(51) 

LWM 

67.3(93) 

66.08(64) 

Normalized Residual 

64.50(50) 

64.33(48) 

Rajeval Technique 

65.57(45) 

65.07(46) 

Best Representation (BR) 

66.2(89) 

66.2(77) 

Bootstrap (BS) 

66.6(15) 

66.0.(14) 

Mandel-Paule 

66.7(92) 

65.3(62) 
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